Atmospheric neutrinos in ice and measurement of neutrino oscillation parameters 
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The main goal of the IceCube Deep Core Array is to search for neutrinos of astrophysical origins. 
Atmospheric neutrinos are commonly considered as a background for these searches. We show that 
the very high statistics atmospheric neutrino data can be used to obtain precise measurements of 
the main oscillation parameters. 



PACS numbers: 14.60.Pq 

I. INTRODUCTION 

A large number of experiments of different types have 
provided strong evidence for neutrino oscillations and 
thus for physics beyond the Standard Model (see Ref. [l[ 
and references therein). 

Cosmic ray interactions in the atmosphere give a nat- 
ural source of neutrinos. These atmospheric neutrinos in 
the GeV range have been used by the Super-Kamiokande 
(SK) detector to provide evidence for neutrino oscilla- 
tions |2j. The large size of neutrino telescopes such as 
AMANDA, IceCube and KM3NeT makes possible the de- 
tection of a large number of atmospheric neutrino events 
with a higher energy threshold, ~ 100 GeV, even though 
the neutrino flux decreases rapidly with energy (~ E~ 3 ). 
These observatories were built to detect neutrinos from 
astrophysical sources, or from the annihilation of Weakly 
Interacting Massive Particles (WIMPs) [3|. Because of 
the high threshold energies, neutrino oscillation effects 
in these ice/ water Cherenkov detectors are small. 

Recently, a low energy extension of the IceCube de- 
tector, the IceCube Deep Core array (ICDC) has been 
proposed and deployed [J] . It consists of six densely in- 
strumented strings (7m spacing between optical modules) 
located in the deep center region of the IceCube detec- 
tor plus the seven nearest standard IceCube strings. Its 
goal is to significantly improve the atmospheric muon 
rejection and to extend the IceCube neutrino detection 
capabilities in the low energy domain, down to muon 
or cascade energies as low as 5 GeV. The instrumented 
volume is 15 Mton. Such a low threshold array bubried 
deep inside IceCube will open up a new energy window on 
the universe. It will search for neutrinos from sources in 
the Southern hemisphere, in particular, from the galactic 
center region, as well as for neutrinos from WIMP annihi- 
lation, as originally motivated. In [5J we have proposed 
neutrino oscillation physics as a further motivation for 
building such an array. In particular, we have analyzed 
the sensitivity of ICDC to the neutrino mass hierarchy. 
ICDC can detect tens of thousands of atmospheric neu- 
trino events per year, orders of magnitude beyond the 
present data sample, providing rich opportunities for de- 
tailed oscillation studies. In |6j we have also analyzed the 
cascade signal in ICDC, showing that this can be used to 



get a high rate of v T interactions. 

Here we concentrate on the muon neutrino disappear- 
ance analysis in ICDC, exploring the precision that can 
be reached in the measurement of the main neutrino os- 
cillation parameters, Am 2 ^ and 6*23. 

Section [TT] reviews briefly our present knowledge of the 
neutrino oscillation parameters as well as their expected 
uncertainties from near future facilities. We then proceed 
in section IIIII to describe the analysis presented here and 
the results obtained from it. We close with an outlook in 
section flV] 

II. NEUTRINO OSCILLATIONS 

Neutrino data from solar, atmospheric, reactor and 
accelerator experiments is well understood in terms of 
three-flavor neutrino oscillations. Two Am 2 values and 
two (large) mixing angles are well determined, while the 
third mixing angle is limited to be very small. The CP- 
violating phase (8) is completely unconstrained. In addi- 
tion, the sign of Am 2 ^ is also unknown. 

The best fit oscillation parameter values obtained from 
present data are Q: 

I Am^ I = 2.56 x 10" 3 eV 2 

Am^ = 7.6 x 10" 5 eV 2 

sin 2 2023 = 0.99 

tan 2 6»i2 = 0.47 (1) 

and sin 2 26»i 3 < 0.15 for Am§ x = 2.5 x 10~ 3 eV 2 . No- 
tice that an extra unknown in the neutrino oscillation 
scenario is the octant in which #23 lies, if sin 2 26*23 7^ 1- 
This has been dubbed in the literature as "the 6*23 octant 
ambiguity". At the three sigma level lAmjjjJ can vary 
between 2.0 — 2.83 x 10~ 3 eV 2 , while 6*23 varies between 
35.5° -53.5°. 

In the near future, long baseline experiments like MI- 
NOS Q and T2K [§| will improve the current precision 
on Am^ and possibly discover a non-zero value of 6*13, 
if this is close to the present upper limit. In a few years, 
reactor experiments like DoubleChooz [l(|, RENO [ll[ 
and Daya Bay [l2j will provide improved sensitivity to 
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013- This information can be used as input in our analy- 
sis, reducing some of the parameter uncertainties. 

In the past, atmospheric neutrinos in the Super- 
Kamiokande detector have shown evidence for neutrino 
oscillations and the first measurements of Am^l and 
sin 2 2#23 [2| j providing compelling evidence for neutrino 
oscillations versus more exotic phenomena [l3j . While 
facing more systematics than accelerator/reactor exper- 
iments due to the uncertainties in the (natural source) 
neutrino fluxes, atmospheric neutrinos provide great op- 
portunities for exploring oscillation physics due to the 
large range of energies and pathlengths that they span. 
The ICDC detector will collect a data sample which is a 
few orders of magnitude larger than that of the Super- 
Kamiokande experiment and can also measure energy 
and directional information, such that many of the sys- 
tematic errors associated with unknown normalizations 
of fluxes, cross-sections, etc. can be much better under- 
stood and reduced by using the data itself. The MIN- 
ERVA [14] experiment will also provide important infor- 
mation about cross-sections in the relevant energy range. 

The downward going neutrinos are largely unaffected 
by oscillations, so they can be used for determining the 
atmospheric neutrino flux. 

In our numerical calculations we have taken into ac- 
count full three flavor oscillations. It is however straight- 
forward to see that solar parameters do not play an im- 
portant role in the analysis due to the rather high energy 
threshold of ICDC. 



III. ANALYSIS AND RESULTS 

We concentrate our analysis on the very high statis- 
tics track signal. This is dominated by muon events 
originating from the charged current (CC) interactions 
of muon neutrinos. Additional small contributions come 
from v e — > oscillations for non-zero values of #13, as 
well as z^ t — >• v T — > r — > fi events. At low energies cas- 
cade events from various interactions of all neutrino fla- 
vors might not be distinguished from tracks, so they also 
contribute to the same event rate. 

We investigate the neutrino energy range between 
10 GeV and 100 GeV, assuming bins of 5 GeV width in 
the observable muon energy. Note that the energy of sec- 
ondary muons from CC interaction in the 10 — 100 GeV 
neutrino energy range of interest here is (E^) = 0.52.E„ 
and (Ep) = 0.66i?p, respectively for neutrinos and an- 
tineutrinos [15]. For our numerical calculations, we use 
differential neutrino interaction cross-sections as given by 
[l5l | with CTEQ6 parton distribution functions. 

As previously mentioned, we only include upward- 
going neutrinos in the analysis for determining the neu- 
trino oscillation parameters. The down-going neutrino 
events are not sensitive to oscillations, but could im- 
prove uncertainties in the normalization of the atmo- 
spheric neutrino flux. 

We evaluate the number of muon events from muon 



neutrino interactions as a function of observable energy 
and zenith angle as: 

N i ^(E^,c v ) = 2irN A pt dE^ dc v 

JEi Jcuj 

f°° d<b d<r cc 

4 dE ^dE%^ ^(^)^W^) 

+v -> v . (2) 

The other contributions to the event rate are obtained 
in a similar way, using the corresponding initial fluxes, 
oscillation probabilities and cross-sections. 

In the above equation, A, and A^ are, respectively, the 
bin widths of the i energy bin and j zenith angular bin 
(defined by c„ = cos 6 ze nith)i t is the exposure time, d<^„'s 
are the atmospheric (anti)ncutrino differential spectra 
and da /dE are differential CC (anti)neutrino cross 
sections. 

The oscillation probabilities P Va _^ v . have been ob- 
tained in a full three flavor scenario, including matter 
effects. While the signal is dominated by muon disap- 
pearance, given the very high statistics, sub-dominant 
effects become observable if #13 is within reach of the 
near future reactor and accelerator experiments. 

For the atmospheric (anti)neutrino fluxes, 
d<\> Va j dE v d£l we use the results from Refs. [l6|. The 
atmospheric neutrino fluxes from Refs. [l7j ] give similar 
results. The absolute electron and muon atmospheric 
(anti) neutrino fluxes have errors of 10% — 15% in the en- 
ergy region of interest here [l8[ . Those errors are mostly 
induced by our incomplete understanding of hadron pro- 
duction, although the situation is expected to improve 
with HARP and MIPP data. In addition, a calibration 
of this overall systematic uncertainty can be performed 
by looking at different angular apertures/energy ranges 
where the oscillatory signal is not present, For the 
neutrino-antineutrino flavor ratios 4>v^/4 > y^ an d 4>vJ4>v c 
the uncertainty is reduced to ~ 7% in the energy range 
we explore in the current study [36| . Even smaller 
uncertainties are expected when the muon-to-electron 
flavor ratio (</>„ + <t>v^)/{<l>v e + </>p e ) is considered. We 
will comment on the impact of the atmospheric neutrino 
flux uncertainties on our results below, when we discuss 
systematic uncertainties in our numerical analysis. 

The factor in Eq. [3 accounts for the effective detec- 
tor volume. Unlike in many other detectors that use fixed 
geometric cuts and sharp detection thresholds, neutrino 
telescopes employ an energy-dependent effective volume, 
which depends on the type of observable considered, in 
order to maximize the number of events in their data 
analysis. Preliminary ICDC estimates Q show that this 
effective volume is of the order of 10 Mtons for 10 GeV 
neutrinos, increasing with energy. We present here re- 
sults for our own effective volume estimates, based on 
several different physics-motivated assumptions regard- 
ing the detector. An optimistic configuration consists 
in using the full ~ 15 Mtons detector mass, correspond- 
ing to the 250 m diameter, 350 m height detector. The 
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conservative option that we have adopted is to consider 
only a 35 m radius around each of the 6 densely instru- 
mented strings (at the ICDC depth, the light attenua- 
tion length is 40-45m). This reduces the total mass to 
~ 7 Mtons, but avoids overestimating the observable low 
energy events, which are crucial for the precise determi- 
nation of oscillation parameters, as we will further discuss 
below. 

For a fixed detector configuration, we considered two 



where R fJ ,(E fl ) is the energy-dependent muon range in 
ice. The fully contained events would yield a better en- 
ergy measurement. On the other hand, the atmospheric 
neutrino flux decreases very fast with energy and retain- 
ing muons that exit the detector helps improve the statis- 
tics at higher energies. These two possibilities turn out 
to give equivalent results in the final analysis, as most 
of the information about oscillation parameters is con- 
tained in the low energy data, where all events are fully 
contained without imposing further constraints. 

We always require the neutrino interaction to occur in- 
side the detector volume. ICDC aims to detect the cas- 
cade from the initial neutrino interaction in addition to 
the muon track, which would provide a better measure- 
ment of the total initial neutrino energy. Availability of 
such information would improve our analysis, which only 
relies on the measurement of the muon energy. 

Since no detailed ICDC simulations of the expected 
systematic uncertainties is presently available, we dis- 
cuss results for different assumptions regarding these un- 
certainties. Systematics that affect the detector observ- 
ables include angular reconstruction, muon track length 
reconstruction (directly related to the energy measure- 
ment) , uncertainties introduced by the modeling of light 
propagation and light detection efficiency of the optical 
modules. The uncertainty on the energy and angular 
measurements is related to the estimation of the effective 
area/volume, which is a function of both the muon en- 
ergy and the muon direction. As previously mentioned, 
ICDC does not have a sharp energy threshold and the 
ICDC analysis shows that the detector can trigger muons 
with energies as low I GeV and reach an effective volume 
of 10 Mtons for fOGeV neutrinos. We use a threshold of 
5G eV muon energy in our analysis. Given the detector 
properties, we estimate the energy resolution to be given 
by a 25m track length uncertainty, which translates into 
about 5 GeV muon energy. 

We explore the sensitivity of the atmospheric neutrino 
data to the main oscillation parameters, Am^ and #23 
using a x 2 statistical analysis, which contains both sta- 



different scenarios to define the effective volume: (1) al- 
low only fully contained events, and (2) allow some of 
the muons exit the detector. In both cases we require 
the initial neutrino interaction occurs inside the detector 
volume. Note that ICDC uses the large IceCube detector 
as veto, in order to exclude events that originate outside 
the ICDC region. For the fully contained events, for a 
detector with cylindrical shape of radius r and height h, 
is given by [l9( 



(3) 

I 

tistical and systematic uncertainties. To account for sys- 
tematic uncertainties we have included an overall error 
computed as a fraction f sys of the number of events in 
each bin: f sys ■ N^j^E^, c v ), added in quadrature with 
the statistical error. We will consider f sys values of 5% 
and 10% for our analysis. This treatment of systematic 
uncertainties is more conservative than other more elab- 
orated approaches used in analysis of neutrino data, as, 
for instance, the so-called pull method [13,[2l|. This sys- 
tematic uncertainty can account for the issues discussed 
above, including uncertainties in fluxes, cross-sections 
and total detector volume. Given that many of the un- 
certainties discussed can be reduced once data becomes 
available and the angular and energy regions not sensitive 
to the oscillation signal are used, a 10% overall system- 
atic uncertainty is likely an overestimate of what can be 
achieved. Other types of uncertainties, concerning an- 
gular resolution, energy resolution and energy thresholds 
are explored in more detail below, as we consider different 
scenarios in our analysis. 

Figure [T] shows the 90% CL contours for the allowed 
parameter values for two degrees of freedom (dof), with 
the atmospheric neutrino data collected in ICDC after a 
10 years exposure. Note that our choice of 2 dof statis- 
tics is rather conservative as explained in the Appendix 
of Ref. [22]. The three plots correspond to three different 
reference values of #23 ■ 45°, 43° and 41°, all within the 1 
sigma presently allowed region. The very large statistics, 
of up to 10,000 events per year, provide very good sen- 
sitivity to oscillation physics and it can be seen that the 
precision will be significantly improved compared with 
the present one. Maximal mixing could be potentially 
excluded for values of 6*23 below 42°. 

The final sensitivity to the parameters of interest is 
however strongly influenced by our knowledge of #13 and 
various types of systematic uncertainties. We explore 
these effects in some detail and show that under reason- 
able assumptions oscillation parameters can be measured 
with much higher precision than present data allows. 

For Figure [T] we have used a reference value of 



VpiE^ff) = 2hr 2 arcsin U 1 - ^^tl^e] ( 1 - MEA | cos , 
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FIG. 1: The left, middle and right panels show the 90% CL contours for 2 dof, for atmospheric mixing angles #23 = 45° , 43° and 
41°, respectively. The different colors account for different energy /angular resolutions. The solid lines allow #13 varying freely 
in its presently allowed range. The dashed lines show the results for a fixed 613 at the assumed true value sin 2 2#i3 = 0.01. The 
dotted lines corresponds to free #13, but with a gaussian prior centered around the true value sin 2 2#i3 = 0.01, with a 1 sigma 
uncertainty of 0.02. A systematic uncertainty of 5% has been included. 



sin 2 2^i3 = 0.01. In order to understand the influence 
of this angle on the final sensitivity, we show results 
for three different assumptions, represented by different 
types of lines. With our present knowledge of #13 the 
results would be those represented by solid lines, which 
correspond to 6> 13 varying freely in its presently allowed 
range. The dashed lines show the results for the unreal- 
istic case of a fixed #13 at the assumed true value. The 
dotted lines corresponds to free #13, but with a gaussian 
prior centered around the true value and with a 1 sigma 
uncertainty of 0.02. This is a very realistic situation, 
which assumes having an independent measurement of 
9i3 from a different experiment, which is expected on the 
relevant time scale. DoubleChooz should reach an un- 
certainty around 0.03 and Daya Bay around 0.01 or even 
smaller, which would improve this analysis even further. 
It can be seen, however, that even with our conserva- 
tive assumption, the sensitivity is very close to the ideal 
situation of fixed #13. 

The different colors on the same plot depict the effects 
of the energy /angular resolution. The blue lines assume 
5 GeV muon energy bins and four different angular bins: 
c u e (-1, -0.9), (-0.9, -0.8), (-0.8, -0.7) and (-0.7, 0). 
The green lines assume that in the lowest energy bins di- 
rectional reconstruction is very hard and a single angular 
bin of c v € (—1,0) is considered. The red lines corre- 
spond to eliminating the two lowest energy bins from the 
data, thus setting the neutrino energy threshold around 
30 GeV. 

It can be seen that poor angular resolution deteriorates 
the results somewhat, but the energy threshold is crucial 
for a good measurement of #23- For astrophysical source 
searches, a 30 GeV energy threshold would be sufficient 
and potentially even desired, in order to reduce the atmo- 
spheric neutrino background. The detector setup allows 
however for a much lower energy threshold, with an esti- 
mated effective volume of the order of 10 MTon at neu- 



trino energies of 10 GeV. The oscillation effects become 
very small above 40 GeV so it is important to analyze the 
data at the lowest energies available in order to obtain 
the maximum amaount of information about neutrino 
oscillation parameters. It is also very important for the 
measurement of atmospheric neutrino parameters to have 
events with energies both above and below the oscillation 
peak and have information about the shape of the energy 
spectrum, which is most sensitive to the oscillation pa- 
rameters. For bins with energy above the peak, the al- 
lowed regions in the (sin# 2 3j ^ m ii) plane curve upwards, 
while for energies smaller than the peak they curve down- 
wards. The combination of the two regions has a small 
intersection and provides a very good measurement. For 
higher energy threshold, for most pathlcngth, only the 
events above the peak would be observed and the sensi- 
tivity to #23 is greatly degraded. 

In addition to the effects of #13, energy and angular 
resolution effects, we have also added a 5% systematic 
error in the above results to account for flux/cross-section 
and other uncertainties, as previously discussed. This 
uncertainty can be greatly reduced by using the data 
itself to fit for a flux parameter and by including down- 
going and higher energy events. 

Figure [2] shows the equivalent to Fig.Q]but with a refer- 
ence value of sin 2 2#i 3 = 0.08. For the larger value of #13 
illustrated in Fig. [5J the contours become slightly asym- 
metric, as expected. This asymmetry in #23 could in prin- 
ciple be used to get sensitivity to the octant. Our analy- 
sis shows that it may be possible to determine the octant 
for small #23) large #13 and favorable assumptions about 
energy/angular resolution and systematics errors. For 
sin 2 2#i3 = 0.08 there is sensitivity to the octant at the 
90% CL if #23 is smaller than 41°, while for sin 2 29 13 = 0.1 
it might be possible to have octant sensitivity for 823 
smaller than 43°. These values are still in the 1 sigma 
presently allowed regions for the atmospheric mixing an- 
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gle #23. The octant degeneracy [23[ is very hard to resolve 
and future neutrino facilities such as neutrino factories 
and/or superbeams [2J-|28|, the combination of super- 
beams with atmospheric neutrino data [29], [30| , or future 
iron calorimeter atmospheric neutrino detectors [3ll | have 
been shown to provide powerful setups to resolve it. High 
statistics atmospheric neutrino data in ICDC could be a 
first step in addressing this degeneracy, if 6*13 turns out 
to be large and 6*23 non-maximal. 

It is important to note that a three flavor analysis con- 
sidering the phyical quantity sin #23 can yield much bet- 
ter sensitivity than a two flavor one in terms of sin 2 26*23. 
In Ref. 32j it has been shown that going from sin 2 26*23 
to sin #23 one can introduce errors of 10 — 20% for 6*23 
values close to 7r/4. 

All the above plots assume a reference value of Am^ = 
2.5 x 10~ 3 eV 2 . Changing the reference value of Am^ es- 
sentially moves the allowed regions up and down, without 
a significant change in the size and shape of the contours. 

Since the track signal is dominated by the muon neu- 
trino disappearance channel, a non-zero CP-violating 
phase has a very small effect on the parameter sensitiv- 
ity. For the largest values of 613, the effect of the phase 
would be to slightly enlarge the Amfi allowed range to- 
wards lower values. 

Considering smaller values for #23 within the 1, 2 or 
3 sigma presently allowed regions makes a better mea- 
surement even easier, allowing for both the exclusion of 
maximal mixing and, for relatively large values of 613 , for 
octant discrimination. We have not shown values of #23 
larger than 45°, as the results for 6*23 = ir/4 + 5 look like 
"mirror images" of those for #23 = 7r/4 — 5. 

Figures [3] and [4] depict the effects of systematic un- 
certainties, which are often a limiting factor in disap- 
pearance experiments. The blue contours assume no sys- 
tematic error, green correspond to 5% and red to 10% 
systematic errors, for nine energy bins and four angu- 
lar bins. In the former two plots we have assumed free 
#13 with a gaussian prior centered around the true value 
sin 2 26*13 = 0-01 and sin 2 26*13 = 0.08 respectively, with a 



1 sigma uncertainty of 0.02. 

Notice that, while systematic errors from flux/cross- 
section uncertainties, etc. affect the parameter measure- 
ment at some level, they still allow for a fairly precise de- 
termination of Am^ and #23. As previously discussed, 
many of the systematics will improve with data, since 
normalization factors can be included as free parameters 
in the fit and determined with good precision, especially 
using the additional energy and angular bins where oscil- 
lations are not important. Information from near future 
experiments can also be used to reduce some of the un- 
certainties. 

Good muon track reconstruction would be extremely 
useful, but relatively precise measurements can be 
achieved even with reduced angular resolution. The most 
important factor in the determination of the main os- 
cillation parameters is the energy threshold, as most of 
the sensitivity comes from neutrino energies between 10 
and 40 GeV. With the 10 GeV neutrino energy thresh- 
old expected to be reached with the 6 string detector 
we have considered here and with realistic assumptions 
about other parameters and uncertainties, the parameter 
determination will be very good. Two additional strings 
will be deployed this year in the center of the detector. 
These will likely improve the quality of the directional 
reconstruction and the energy threshold and resolution 
and will thus lead to even better results. 



IV. OUTLOOK 

The IceCube detector and its Deep Core array provide 
an excellent opportunity for studies of atmospheric neu- 
trinos. Being the largest existing neutrino detector, it 
will accumulate a huge number of atmospheric neutrino 
events over an enormous energy range, thus allowing for 
detailed studies of oscillation physics, Earth density, at- 
mospheric neutrino fluxes and new physics [33| . In order 
to extract all this information it is necessary to use en- 
ergy and angular distribution information, as well as fla- 
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FIG. 3: The left, middle and right panels show the 90% CL contours for 2 dof assuming #13 = 0.01 and #23 = 45°, 43° and and 
41°, respectively. The different colors account for different values of the systematic uncertainties, see the text for details. We 
have assumed free 813 with a gaussian prior centered around the true value sin 2 2#i3 = 0.01, with a 1 sigma uncertainty of 0.02. 
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vor composition, all possible to obtain with the IceCube 
detector. 

Qualitatively, there are three main energy intervals and 
three main angular regions which are sensitive to different 
types of physics. 

At very high energies, above 10 TeV, neutrino inter- 
action cross-sections become high enough that neutrinos 
going through the Earth start being attenuated 34] . This 
effect is sensitive to neutrino interaction cross-sections 
and to the density profile of the Earth. 

The "intermediate" energy region, between 50 GeV 
and 1 TeV can provide good information about the at- 
mospheric neutrino flux, which can be used to improve 
the uncertainties in the simulated atmospheric neutrino 
fluxes [351 . 

In our paper we concentrated on the "low" energy re- 
gion, below about 40 GeV, where neutrino oscillation ef- 
fects can be significant. Although the IceCube detector 
has higher energy threshold, its low energy extension, the 
Ice Cube Deep Core array (ICDC), extends the IceCube 
neutrino detection capabilities in the low energy domain, 
down to muon or cascade energies as low as 5 GeV. 



This data can be used for precise measurements of os- 
cillation parameters like Am^ and O23 and for resolving 
some of the neutrino oscillation parameter degeneracies. 

A very large effort is directed toward the optimization 
and construction of a next generation of experiments that 
will address precision measurements and potential dis- 
covery of new phenomena in neutrino oscillations. While 
they will achieve extremely good sensitivities, these ex- 
periments face long construction and data-taking times. 
In the meantime, ICDC, which is already taking data, 
will acquire a high statistics atmospheric neutrino sam- 
ple which can be very useful in extracting information 
about neutrino oscillation physics. Some of this infor- 
mation could be used for further optimization of the fu- 
ture accelerator experiments. Long baseline experiments 
can achieve very high precision, but have a fixed baseline 
and a limitted energy coverage. Atmospheric neutrinos 
in ICDC cover a much more extended range of energies 
and many baselines and thus provide potentially new, 
complementary information. Combining data from the 
two types of experiments would check the consistency of 
allowed neutrino oscillation parameter space over a large 
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range of energies and propagation distances and can be 
very important for solving parameter degeneracies. If 
any new physics is present in the neutrino sector, its ef- 
fects could be larger at high energies and long distances 
and would lead to additional effects in the atmospheric 
neutrino data. 

In summary, ICDC offers a unique window toward a 
better understanding of neutrino properties due to its 
very high atmospheric neutrino statistics. Careful stud- 
ies of the expected atmospheric neutrino oscillation sig- 
nals in ICDC, as the one carried out here, are extremely 
important, since atmospheric neutrinos constitute an ir- 
reducible background to astrophysical neutrino searches 
and can offer information that is complementary to that 
obtained from other neutrino experiments. 
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